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ABSTRACT 

We  analyzed  the  full  Stokes  spectra  using  simultaneous  measurements  of  the  photospheric  (Fe  I  630.15 
and  630.25  nm)  and  chromospheric  (Mg  I  hi  517.27  nm)  lines.  The  data  were  obtained  with  the  HAO/NSO 
Advanced  Stokes  Polarimeter,  about  a  near  disc  center  sunspot  region,  NOAA  AR  9661.  We  compare  the 
characteristics  of  Stokes  profiles  in  terms  of  Doppler  shifts  and  asymmetries  among  the  three  spectral  lines, 
which  helps  us  to  better  understand  the  chromospheric  lines  and  the  magnetic  and  flow  fields  in  different 
magnetic  regions.  The  main  results  are:  (1)  For  penumbral  area  observed  by  the  photospheric  Fe  I  lines, 
Doppler  velocities  derived  from  Stokes  I  (vj)  are  very  close  to  those  derived  from  linear  polarization  profiles 
{vip)  but  significantly  different  from  those  derived  from  Stokes  V  profiles  (vzc),  which  provides  direct  and  strong 
evidence  that  the  penumbral  Evershed  flows  are  magnetized  and  mainly  carried  by  the  horizontal  magnetic 
component.  (2)  The  rudimentary  inverse  Evershed  effect  observed  by  the  Mg  I  bi  line  provides  a  qualitative 
evidence  on  its  formation  height  that  is  around  or  just  above  the  temperature  minimum  region.  (3)  vzc  and  vip 
in  penumbrae  and  vzc  in  pores  generally  approach  their  u,  observed  by  the  chromospheric  Mg  I  line,  which 
is  not  the  case  for  the  photospheric  Fe  I  lines.  (4)  Outer  penumbrae  and  pores  show  similar  behavior  of  the 
Stokes  V  asymmetries  that  tend  to  change  from  positive  values  in  the  photosphere  (Fe  I  lines)  to  negative 
values  in  the  low  chromosphere  (Mg  I  line).  (5)  The  Stokes  V  profiles  in  plage  regions  are  highly  asymmetric 
in  the  photosphere  and  more  symmetric  in  the  low  chromosphere.  (6)  Strong  red  shifts  and  large  asymmetries 
are  found  around  the  magnetic  polarity  inversion  line  within  the  common  penumbra  of  the  S  spot.  We  offer 
explanations  or  speculations  to  the  observed  discrepancies  between  the  photospheric  and  chromospheric  lines 
in  terms  of  the  three-dimensional  structure  of  the  magnetic  and  velocity  fields.  This  study  thus  emphasizes  the 
importance  of  spectro-polarimetry  using  chromospheric  lines. 

Subject  headings:  Sun:  activity  —  Sun:  atmospheric  motions  —  Sun:  chromosphere  —  Sun:  magnetic  fields 
—  Sun:  photosphere  —  sunspots  —  line:  profiles  —  techniques:  polarimetric 


1.  INTRODUCTION 

Solar  photospheric  magnetic  fields  and  their  associated 
dynamics  have  been  extensively  studied.  However  an 
understanding  of  their  vertical  variation  from  the  pho¬ 
tosphere  though  the  chromosphere  to  the  corona  needs 
better  comprehension.  Direct  measurement  and  diagnosis 
of  magnetic  and  flow  fields  in  higher  layers,  especially  in 
the  near  force-free  chromosphere  can  play  an  important 
role  in  fully  understanding  their  three-dimensional  (3D) 
structures.  Nevertheless,  such  efforts  were  relatively  rare  due 
to  the  paucity  of  chromospheric  spectral  lines  with  suitable 
Zeeman-split  sensitivity  dDalgarno  &  Lavzerl  i!987|)  and 
the  complicated  dynamics  and  topology  of  the  magnetized 
plasma  in  this  particular  layer.  Taking  advantage  of  modern 
observational  techniques  and  recognizing  the  critical  need  to 
explore  chromospheric  magnetic  fields,  increasing  attention 
has  recently  been  paid  to  the  spectro-polarimetry  of  the  chro 


mosphere  (e.g„  iSocas-Navarro  et  alJ  120001:  IZhang  &  Zhang 


2000;  Trujillo  Bueno  et  al.  2002t  Lopez  Ariste  &  Casini 


2002;  [Solanki  et  al.  2003;  Leka  &  Metcalf!  120031;  lLagg  et  aT 


2004:  Balasubramaniam  et  alJ  120041),  which  reveal  the 


height-dependent  variation  of  different  types  of  magnetic 
features.  Those  analyses  were  made  using  magnetically 
sensitive  spectral  lines  formed  in  the  chromospheric  levels. 


such  as  Ca  II  849.8  and  854.2  nm.  He  I  1083.0  nm  and  D3 
line,  Na  I  Di  589.6  nm,  Ha  and  H  i  lines.  In  particular,  the 
Mg  1 1>2  517.27  nm  line  is  favorable  for  spectro-polarimetry 
in  the  low  chromosphere  because  its  core  forms  in  a  narrow 
region  near  or  right  above  the  temperature  minimum  region 
(TMR)  with  a  relatively  large  Lande  gc//  factor  of  1.75 
(e.g„  iBriand  &  Solankil  fl995h.  Theoretical  calculation 
( Altrock  &  Canfieldl  119741;  iLitesetal]  Il988t  iMauasetahl 


and  observational  analysis  (Briand  &  Solankil  199 


Briand  &  Martinez  Pilled  l200ll  iGosain  &  Choudhary  120031 


have  been  carried  out  on  Mg  I  tu  Stokes  polarization  spectra, 
which  corroborate  that  it  is  well  suited  for  probing  the 
magnetic  and  flow  fields  in  the  low  chromosphere. 

Doppler  shifts  of  Stokes  I  and  Q,U,V  signals  can  be  used 
to  investigate  mass  flows  in  and  around  magnetic  elements 
(ISolankil Il986h  ISolanki  &  Pahlke  1988).  In  contrast  to  the 
Stokes  I  that  represents  the  total  intensity  of  light,  the  Stokes 
Q,U  (the  linearly  polarized  component)  and  V  (the  circularly 
polarized  component)  have  contribution  only  from  the  magne¬ 
tized  plasma.  Thus  Doppler  shifts  of  the  Stokes  I  and  Q ,  U,  V 
profiles  reflect  average  line-of-sight  (LOS)  velocities  of  plas¬ 
mas  in  the  whole  resolution  element  and  those  within  mag¬ 
netic  regions,  respectively  (e.gJSolankill  1986b. 

Asymmetries  of  Stokes  profiles  are  a  powerful  diagnos- 
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DENG  ET  AL. 


tics  of  the  height  and  spatial  properties  of  magnetic  and 
flow  fields  (e.g..  iGrigorev  &  Katzl  119751).  Under  a  restric¬ 
tive  Milne-Eddington  atmosphere  where  magnetohydrody¬ 
namic  (MHD)  conditions  (e.g.,  plasma  velocity,  magnetic 
field  vector,  and  temperature)  are  constant  with  depth  and 
the  source  function  varies  linearly  with  optical  depth  jjUnno| 

ll956t(Landi  Deglinnocenti  &  Landi  Deglinnocentil  19851).  the 

Zeeman-split  Stokes  I,Q,U  profiles  are  symmetric,  while  V 
is  antisymmetric  about  the  line  center.  Most  of  the  observed 
Stokes  spectra,  however,  deviate  from  the  ideal  forms  and 
exhibit  notable  difference  between  the  blue  and  red  wings 
in  amplitude  and/or  area.  Therefore  the  Stokes  asymme¬ 
try  is  able  to  reveal  the  inhomogeneous  nature  of  the  so 
lar  atmosphere  (e.g..  Illling  et  al.l  1974a, b).  Direct  observa 


tions  (e.g.,  Balasubramaniam  et  al.  1997;  Sigwarth  2001)  and 


1988 


1993 


modeling  of  synthetic  Stokes  profiles  (e.g. ,  [Solanki  &  Pahlke 
Sanchez  Almeid^&LitesL^99^JSrfanki_&Montavoij 
Sanchez  Almeida  et  al.  1996;  Grossmann-Doerth  et  al 


2000;  Steinerll2000t)  have  found  possible  physical  conditions 
for  producing  asymmetric  Stokes  spectra,  such  as  gradients 
of  flow  velocity  and/or  magnetic  field  vector  along  the  LOS 
direction  and  nonuniform  atmosphere  that  contains  two  or 
more  distinct  magnetic  and/or  flow  components_in  one  reso¬ 
lution  element.  A  more  comprehensive  study  bvlLopez  Aristel 
(120021)  based  on  the  analytical  solution  for  the  generalized 
transfer  equation  explicitly  demonstrates  that  the  velocity  gra¬ 
dients  are  the  necessary  and  sufficient  condition  for  Stokes 
profile  asymmetries.  The  addition  of  other  effects  or  gradi¬ 
ents  of  other  MHD  parameters  (e.g.,  magnetic  field  and  tem¬ 
perature)  could  enhance  the  asymmetries  and  give  different 
characteristics  of  asymmetries.  For  example,  strong  Stokes  I 
and  V  asymmetries  are  usually  found  around  magnetic  flux 
concentrations  in  the  photosphere,  which  is  attributed  to  the 
presence  of  the  “magnetic  canopy”  (|Grossmann-Doerth  et  all 


1 9SS.fEMlSolankilil  989ilLeka  &  Steineil2001l).  Strong  gra¬ 

dients  in  both  magnetic  and  flow  fields  are  naturally  generated 
because  the  “magnetic  canopy”  spatially  separates  the  field- 
free  convective  downdrafts  from  the  more  static  flux  tube  that 
fans  out  above  the  photosphere. 

The  aforementioned  analyses  of  Stokes  profiles  predomi¬ 
nantly  use  photospheric  spectral  lines.  In  retrospect,  only  few 
such  attempts  using  chromospheric  Stokes  spectra  have  been 
reported.  As  examples. iBriand  &  Solanki  (i  19981)  inferred  fast 
down  flows  concentrated  in  narrow  lanes  surrounding  the  flux 
tube  below  the  “magnetic  canopy”  by  using  Stokes  V  asym- 
metries_of  Mg  I  bi  line  in  additional  to  photospheric  lines. 
iGosain  &  Choudharvl  (2003')  compared  the  Stokes  V  asym¬ 
metries  between  photospheric  Fe  I  lines  and  chromospheric 
Mg  I  lines  and  found  differences  in  their  spatial  distribution 
within  a  simple  sunspot.  Additional  investigation  of  chromo¬ 
spheric  Stokes  spectra  in  conjunction  with  photospheric  mea¬ 
surements  would  improve  our  understanding  of  3D  magnetic 
topology  and  plasma  flows. 

In  this  paper,  we  study  Doppler  shifts  and  asymmetries 
of  Stokes  profiles  of  three  spectral  lines  simultaneously  ob¬ 
served  in  an  active  region  close  to  the  disk  center.  The  three 
spectral  lines,  Fe  I  630.25  nm  (ge//  =  2.5),  Fe  I  630.15  nm 
(geff  =  1  -67),  and  Mg  I  bo,  are  formed  at  partially  overlap¬ 
ping  atmospheric  layers  spanning  from  the  photosphere  to  the 
low  chromosphere.  Even  though  both  are  regarded  as  photo¬ 
spheric  lines,  the  Fe  I  630.15  nm  line  often  forms  above  the 
630.25  nm  line  (e.g..lKhomenko  &  Colladosll2007l).  Combin¬ 
ing  data  from  multiple  atmospheric  heights  using  these  spec¬ 
tral  lines  thus  enables  us  to  shed  new  light  on  the  properties  of 
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Fig.  1 . —  ( a )  The  Fe  I  630. 15  nm  continuum  image  of  the  NOAA  AR  9661 
on  2001  October  17.  White  and  black  contours  outline  the  positive  and  neg¬ 
ative  polarities,  respectively.  ( b )  The  corresponding  longitudinal  magnetic 
field  mimicked  by  the  peak  value  of  Stokes  V  blue  lobe. 


magnetic  and  flow  structures  associated  with  different  mag¬ 
netic  features.  The  plan  of  this  paper  is  as  follows:  in  §§|2] 
and  [TJ  we  introduce  the  observations  and  present  the  proce¬ 
dure  for  data  reduction.  In  §  [T]  the  results  of  comparison  of 
profile  shifts  and  asymmetries  among  the  three  spectral  lines 
are  described  and  discussed.  We  summarize  our  major  find¬ 
ings  and  present  some  prospects  for  future  studies  in  §0 

2.  OBSERVATIONS 

We  obtained  the  spectropolarimetric  data  sets  using  the 
76  cm  Dunn  Solar  Telescope  and  its  associated  Adaptive  Op¬ 
tics  system  (tRimmelei  [2004).  along  with  the  High  Altitude 
Observatory /National  Solar  Observatory  Advanced  Stokes 
Polarimeter  (HAO/NSO  ASP.lElmore  et  al.ll  1 9921:  lLitesHT996l: 
Skumanich  et  all 1 19971).  The  observations  were  obtained  on 
2001  October  17  under  excellent  seeing  conditions.  The 
ASP  is  a  slit-grating  spectropolarimeter.  The  1.6'  x  0.6"  slit 
was  oriented  in  the  solar  north-south  direction  and  scanned 
in  the  east-west  direction.  Each  whole  scan  took  280  steps 
with  a  step  size  of  0.6",  which  results  in  a  field  of  view 
(FOV)  of  1.6'  x  2.8'  that  covers  the  near  disk  center  (N14°, 
W6°,  /;  =cos0  =  0.96)  active  region  NOAA  9661  as  shown 
in  Figure  [l]  It  contains  a  S  sunspot  at  left,  an  a  sunspot  at 
right,  as  well  as  several  pores  and  plages  lying  in  between. 
The  obtained  spatial  sampling  is  about  0.42"  pixel-1  in  the 
north-south  direction  and  0.6"  pixel-1  in  the  east-west  direc¬ 
tion.  Full  Stokes  I,Q,U,V  spectra  were  taken  simultaneously 
in  two  spectral  bands  centered  at  630.2  nm  and  517.27  nm 
by  two  ASP  cameras.  The  dispersions  for  the  two  spectral 
bands  are  1.27  and  1.02  pm  pixel-1,  respectively.  Each  whole 
scan  took  about  25  minutes.  A  total  of  five  whole  scans 
were  continuously  taken  within  2.5  hours,  during  which  no 
significant  evolution  of  the  magnetic  structures  was  found. 
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Other  information  about  the  observation  run  can  be  found  in 
IChoudharv  &  Balasubramaniaml  (120071). 


3.  DATA  REDUCTION 

We  applied  the  standard  ASP  calibration  procedures 
dLites  et  alll991brSkumanich  et  al.l  19971)  to  the  data  sets.  The 
calibrated  Stokes  I,Q,U  ,V  spectra  were  normalized  to  the 
quiet  Sun  continuum  intensity  Ic  that  was  obtained  by  fitting 
the  Kitt  Peak  FTS  atlas  to  the  observed  Stokes  I  profiles.  The 
data  sets  from  the  five  whole  scans  were  then  integrated  to  in¬ 
crease  the  signal  to  noise  ratio  (S/N).  Furthermore,  we  calcu¬ 
lated  the  standard  deviation  a  in  near  continuum  wavelength 
range  for  each  integrated  Stokes  /,  Q.  U,  V  profile  to  represent 
the  profile  noise  level.  Only  profiles  with  signal  amplitude 
greater  than  la  (an  empirical  S/N  threshold  that  performs  well 
in  our  data  analysis)  were  used  for  this  study. 

We  extracted  the  parameters  for  analyzing  the  shift  and 
asymmetry  of  Stokes  profiles  using  the  following  procedures. 

(1)  Stokes  I  LOS  velocities  Vj.  We  employed  a  Fourier  phase 
method  (see  ISchmidt  et  al.l  f 1999’)  to  determine  the  line  core 
wavelength.  We  then  converted  the  wavelength  shifts  to  ve¬ 
locities  (i/j)  by  using  the  Doppler  formula.  As  a  frame  of  refer¬ 
ence,  we  use  the  average  velocity  of  a  small  area  in  the  nearly 
motionless  umbra. 

(2)  Stokes  V  LOS  velocities  vzc.  Normalized  Stokes  V  pro¬ 
files  were  smoothed  to  remove  local  noise.  They  were  then 
divided  into  two  cases:  normal  profiles  that  have  two  opposite 
lobes  and  one  zero-crossing  (ZC)  point  in  between,  and  abnor¬ 
mal  profiles  that  have  no  or  multiple  ZCs  within  a  wavelength 
range  corresponding  to  ±5  kms-1  from  Stokes  I  line  center. 
Using  this  method,  there  is  a  small  chance  that  we  may  im¬ 
properly  classify  normal  profiles  into  an  abnormal  class.  For 
example,  strong  magneto-optical  effects  (see  Fig.  [3k  and  text 
in  §  14.  Ill  produce  Stokes  V  profiles  with  3  ZCs,  although  the 
profile  is  essentially  normal.  Such  improperly  classified  cases 
only  account  for  a  small  portion  (<  0.3%)  in  the  data  and  do 
not  affect  the  analysis.  For  normal  profiles,  the  ZC  shifts  are 
converted  to  LOS  velocities  (yzc)  in  a  same  manner  as  //,.  The 
abnormal  profiles  will  be  discussed  in  §14.  II 

(3)  Linear  polarization  (LP)  LOS  velocities  uip.  The  com¬ 
bination  ( Q2  +  U2)1/2  is_a  measure  of  the  overall  LP  magnitude 
(e.g..lLeka  &  Steinerll2001).  The  LP  profiles  were  smoothed 
to  remove  local  noise.  It  is  difficult  to  measure  the  position  of 
the  central  n  component  as  it  is  usually  small  and  sometimes 
even  vanishes  (see  Figure  [7]).  Instead  we  measure  the  peak 
positions  of  the  two  a  components  and  use  their  mid-point 
as  the  position  of  the  LP  profile.  We  then  converted  the  LP 
profile  shifts  to  LOS  velocities  (vip)  in  a  same  manner  as  v,. 

(4)  Amplitude  asymmetry  8a  and  area  asymmetry  5A.  The 
peak  amplitudes  of  the  blue  (at,)  and  the  red  (a,-)  lobes  of  nor¬ 
mal  Stokes  V  profiles  and  those  of  the  two  a  components  of 
LP  profiles  are  determined.  The  areas  of  blue  (At,)  and  red 
( Ar )  lobes  are  obtained  by  integrating  over  the  same  wave¬ 
length  range  (0.02  nm)  on  both  sides  of  the  ZC  wavelength  for 
Stokes  V  profiles  or  the  mid-point  of  the  two  a  components 
for  LP  profiles.  The  integrating  range  not  only  covers  the  ma¬ 
jority  of  the  Stokes  V  or  LP  signal  corresponding  to  Stokes 
I  Doppler  core,  but  also  rules  out  the  blends  from  the  neigh- 
boring  telluric  lines.  FollowinglSolanki  &  Stenflol(ll984l).  the 
amplitude  asymmetry  8a  and  area  asymmetry  8A  are  derived 
using  the  following  formulae: 


s  kl-kl 

0a  =  - - ; - : - r 


(i) 


\Ab\  +  \Ar\  ■ 


(2) 


4.  RESULTS  AND  DISCUSSIONS 

Figure  [2]  shows  the  Stokes  /,  V  and  LP  images  constructed 
by  line  core  intensity  and  integrated  Stokes  V  and  LP  sig¬ 
nals,  respectively.  Images  of  the  three  spectral  lines  are  ar¬ 
ranged  following  the  order  of  their  formation  heights  (increas¬ 
ing  height  from  bottom  to  top).  The  orange  points  are  the  lo¬ 
cations  of  noisy  profiles  with  S/N  <  7.  The  LP  spectra  have 
more  noisy  profiles  than  Stokes  V,  because  the  Zeeman  sig¬ 
nals  for  LP  are  essentially  lower  and  near  the  disk  center  the 
transverse  magnetic  fields  are  weaker  than  the  longitudinal 
fields  for  most  of  the  areas  in  the  low  atmosphere.  Compared 
to  the  photospheric  Fe  I  spectra,  the  chromospheric  Mg  I  po- 
larimetric  (Stokes  V  and  LP)  spectra  generally  have  lower  S/N 
thus  more  useless  profiles.  The  reasons  are  mainly  due  to  the 
larger  line  width  and  the  lower  magnetic  field  strengths  in  the 
chromosphere.  In  the  following,  we  first  address  the  spatial 
distribution  of  abnormal  Stokes  V  profiles,  which  are  impor¬ 
tant  for  evaluating  the  complexity  of  magnetic  and  flow  fields 
(e.g..  lMickevlll985l:  ISigwarthll200lL  and  references  therein). 
We  then  present  the  results  of  Doppler  shifts  and  asymmet¬ 
ric  properties  of  the  Stokes  profiles,  where  all  the  noisy  and 
abnormal  profiles  are  excluded  from  analysis. 

4. 1 .  Abnormal  Stokes  V  Profiles 

In  the  middle  column  of  Figure  [2j  locations  with  abnor¬ 
mal  Stokes  V  profiles  are  marked  using  colored  points.  They 
are  divided  into  three  kinds  according  to  the  number  of  ZC 
within  a  range  of  ±5  kms-1  from  their  Stokes  I  line  cen¬ 
ters  as  shown  in  Figure 0  Profiles  with  two  ZCs  (green)  are 
most  frequently  seen,  which  have  three  lobes  with  the  cen¬ 
tral  one  opposite  to  the  other  two  (Fig.  0?).  They  concen¬ 
trate  in  the  regions  of  magnetic  polarity  inversion  line  (MPIL) 
and  the  outer  penumbral  boundaries.  They  are  mainly  caused 
by  the  presence  of  mixed  polarities  associated  with  differ¬ 
ent  velocities  (e.g..  ISigwarthl  12001).  whereby  the  two  polar¬ 
ities  can  be  present  beside  each  other  (i.e.,  within  one  res¬ 
olution  element)  or  along  the  LOS  over  the  line  formation 
region  (ISanchez  Almeida  &  Litesll992l:lSolanki  &  Montavonl 
1 19931).  In  outer  penumbral  regions,  the  number  of  two  ZCs 
decreases  with  height,  which  implies  that  the  mixed  polarity 
effect  is  weaker  at  higher  atmospheric  levels.  This  is  prob¬ 
ably  due  to  the  fact  that  the  orientation  of  outer  penumbral 
fields  changes  from  horizontal  or  even  downward  to  more  ver¬ 
tical  when  they  extend  from  the  photosphere  to  the  chromo¬ 
sphere  (IChoudharv  &  Balasubramaniaml2007L  and  references 
therein).  However,  in  the  central  part  of  the  MPIL  of  the  8 
spot  there  are  more  locations  with  two  ZCs  in  the  Mg  I  image. 
We  speculate  that  this  may  point  to  a  stronger  mixed  polarity 
effect  around  the  highly  non-potential  MPIL  in  the  chromo¬ 
sphere. 

The  central  part  of  the  <5  spot’s  MPIL  is  also  populated  by 
locations  without  ZC  (red).  Such  locations  may  involve  mag¬ 
netic  components  carrying  high  speed  flows,  which  cause  the 
ZC  to  shift  far  from  the  line  center  (see  Fig.[3j/).  Since  the  ZC 
position  of  Stokes  V  is  susceptible  to  the  presence  of  veloc¬ 
ity  gradients  through  the  line  formation  region  dLopez  Aristel 
120021)  or  distortion  of  line  profiles  by  noise,  it  is  also  possible 
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Fig.  2. —  Stokes  I  (line  core  intensity),  Stokes  V  (integration  over  Stokes  V  blue  lobe)  and  LP  (integration  over  LP  signal)  images  in  the  three  spectral  lines. 
The  sunspot  umbrae,  penumbrae  and  pores,  as  identified  on  the  basis  of  continuum  intensity,  are  outlined  by  the  yellow  contours.  The  orange  points  on  Stokes  /, 
V  and  LP  images  indicate  the  locations  where  the  Stokes  7,  V  and  LP  profile  has  a  S/N  less  than  7,  respectively.  The  positions  with  abnormal  Stokes  V  profiles 
(red,  green,  and  blue  for  zero,  two,  and  three  ZCs,  respectively)  are  also  denoted  on  each  Stokes  V  image  accordingly.  To  enhance  the  visibility,  the  size  of  the 
colored  points  is  larger  than  the  real  size  of  the  corresponding  image  pixel. 


that  the  profiles  at  these  locations  are  highly  asymmetric  or 
distorted. 

There  exists  another  kind  of  profile  with  three  ZCs  (blue). 
They  could  be  caused  by  residual  noise  or  mixed  polarity  ef¬ 
fect.  Other  than  this,  the  three  ZCs  profiles  within  the  umbra 
seen  in  630.25  nm  Stokes  V  image  are  most  likely  due  to  the 
magneto-optical  effect  (see  Fig.0j:)  caused  by  strong  Zeeman 
effect  (e.g.JWest  &  Hagyardll  19831). 

4.2.  Doppler  Shifts  of  Stokes  I,  V  and  LP  Profiles 

In  the  upper  panels  of  Figure®  we  show  Dopplergrams  (z/,-, 
vzc  and  vip)  for  each  spectral  line.  They  are  constructed  using 
the  LOS  velocities  derived  from  the  shifts  of  Stokes  I,  V  and 
LP  profiles  at  each  spatial  point,  respectively.  We  note  the 
followings:  (1)  For  the  two  photospheric  lines,  the  spatial  dis¬ 
tribution  of  penumbral  Evershed  flow  in  is,  Dopplergrams  is 
similar  to  that  in  isip  but  significantly  different  from  vzc  Dopp¬ 
lergrams.  (2)  The  variation  of  Evershed  flow  with  height  can 
be  seen  by  the  three  spectral  lines.  Interestingly,  the  Ever¬ 
shed  flow  reverses  its  sign  from  the  Fe  I  to  Mg  I  Doppler¬ 
grams  in  some  places  as  illustrated  by  a  center-side  penumbra 
(white  box),  which  is  known  as  the  inverse  Evershed  effect 
observed  in  the  chromosphere  (e.g..lSolankill2003).  While  in 
some  other  places,  the  sign  reversal  has  not  yet  happened  or 
is  rudimentary.  This  provides  a  qualitative  evidence  on  the 
formation  height  of  the  Mg  I  b2,  which  is  at  or  just  above  the 
TMR.  (3)  The  Fe  I  is,  Dopplergrams  show  that  small  pores  and 
plages  lying  in  between  the  two  sunspots  are  generally  sur¬ 


rounded  by  red-shifted  downflows,  some  of  which  do  not  ap¬ 
pear  in  the  Mg  I  is,  Dopplergram  as  pointed  by  white  arrows  on 
the  top  three  panels  of  Figure0]  In  contrast,  the  Fe  I  iszc  Dopp¬ 
lergrams  show  ZC  redshifts  all  over  the  plage  regions.  We 
suspect  that  the  red-shifted  vzc  signal  may  not  represent  real 
down  flows  in  flux  tubes  but  could  result  from  the  combina¬ 
tion  of  the  susceptibleness  of  ZC  position  to  the  large  Stokes 
V  asymmetry  often  observed  by  photospheric  lines  in  plage 
regions  and  the  limited  spectral  resolution  due  to  smoothing 
(see  ISolanki  &  Sfenfloiri986t  iLonez  Aristell2002l) . 

We  quantitatively  compare  the  is,  with  iszc  and  vip  in  differ¬ 
ent  magnetic  regions  for  each  spectral  line  by  making  scatter 
plots  of  is,  vs.  iszc  and  is,  vs.  is/p  shown  in  the  forth  and  fifth 
row  of  Figure  01  respectively.  Since  the  LP  signal  is  only 
strong  enough  in  penumbrae  and  MPIL,  velocities  from  these 
regions  are  plot  in  isj  vs.  is/p  panels.  Each  group  of  data  points 
are  fit  with  a  form  of  iszc  =  kisj  +  b  or  isip  =  kis ,•  +  b,  and  their 
linear  correlation  coefficients  (C.C.)  are  also  obtained.  We 
discuss  the  scatter  plots  for  different  regions  as  follows. 

Pores — For  the  two  photospheric  Fe  I  lines,  \iszc\  <  (i.e., 

k  <  1)  is  observed  with  moderate  correlation  (C.C.=  0.56  and 
0.70,  respectively).  For  the  Mg  I  line,  | iszc |  «  |i/,j  (i.e.,  k  ss  1) 
with  strong  correlation  (C.C.=  0.85)  is  observed.  It  is  well 
established  that  the  magnetic  filling  factor  (MFF)  is  almost  1 
and  the  magnetic  fields  are  nearly  vertical  in  the  pores.  This 
is  even  true  for  pores  in  heights  above  the  photosphere,  which 
explains  that  iszc  show  large  similarity  with  is-,  as  seen  by  the 
Mg  I  line.  While  down  to  the  photosphere  near  the  edge  of 
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Fig.  3. —  The  sample  abnormal  Stokes  V  profiles  defined  in  this  study  with  zero,  two,  and  three  ZCs  within  ±5  kms  1  from  their  Stokes  I  line  centers  indicated 
by  the  vertical  solid  line.  The  vertical  dotted  lines  mark  the  wavelength  range  corresponding  to  ±5  kms-1  from  Stokes  I  line  center. 


pores  the  field  lines  may  be  more  inclined  and  the  MFF  may 
be  less  than  1  in  such  an  observation  with  median  spatial  res¬ 
olution  (e.g.,  Leka  &  Steiner  2001).  Strong  and  narrow  field- 
free  down  flows  with  speed  larger  than  that  within  the  pores 
have  been  found  immediately  around  pores  in  the  photosphere 
(tSankarasubramanian  &  Rimmeleil2003l).  which  gives  rise  to 
|  Vi  |  instead  of  \vzc\.  This  might  be  responsible  for  the  observed 
\vzc\  <  Wi\  in  the  Fe  I  lines. 

Penumbrae — The  penumbral  magnetic  fields  are  known  to 
consist  of  a  more  vertical  component  and  a  more  horizon¬ 
tal  component  in  the  photosphere  (interlocking-comb  struc¬ 
ture,  e.g.,  llchimoto  et  al.l 12007).  For  this  region  close  to  the 
disk  center,  vzc  have  more  contribution  from  the  more  verti¬ 
cal  component  than  from  the  other  component.  On  the  con¬ 
trary,  vip  have  more  contribution  from  the  more  horizontal 
component.  Figure  Q]  shows  that  v,  has  much  better  spatial 
and  magnitude  correlation  with  vip  than  with  vzc  for  the  two 
Fe  I  lines.  Moreover,  the  slopes  of  v,  vs.  vip  (k=  0.87  and 
0.86,  respectively)  are  close  to  1  and  much  larger  than  those 
of  Vj  vs.  vzc  (k=  0.33  and  0.35,  respectively).  In  other  words, 
the  vip  is  very  close  to  v,  and  their  magnitude  is  much  larger 
than  |  vzc  | .  These  results  directly  support  that  the  penumbral 
Evershed  flows  are  not  only  magnetized  but  mainly  carried 
by  the  horizontal  magnetic  component  as  previously  found  by 
other  authors  using  different  observations  and  methods  (see, 
e.g..lRimmele  &  Marinoll2006:  Ichimoto  et  alJl2008l  and  ref¬ 
erences  therein).  While  for  the  Mg  I  line,  Vj  is  observed  sim¬ 
ilar  to  both  of  vip  and  vzc  (i.e.,  k  «  1)  with  small  correla¬ 
tion  coefficients  (C.C.=  0.21  and  0.40,  respectively).  This 
implies  that  rising  up  to  the  TMR  the  penumbral  magnetic 
and  flow  structures  are  different  from  those  in  the  photo¬ 
sphere.  The  large  scatter  of  the  data  points  resulted  from  nois¬ 
ier  and  weaker  polarimetric  signal  of  the  Mg  I  line  hinders  us 
from  drawing  solid  conclusion.  We  thus  speculate  that  the 
interlocking-comb  structure  of  the  penumbral  magnetic  and 
flow  field  might  be  a  shallow  phenomenon  mainly  in  the  pho¬ 
tosphere. 

MPIL — For  a  compact  elongated  region  around  the  MPIL 
and  within  the  common  penumbra  of  the  8  spot,  apart  from 
noise  dominating  the  Fe  I  630.25  Vj  and  v/p  signal,  the  Dopp- 
lergrams  show  strong  redshift  with  a  mean  speed  conspicu¬ 


ously  larger  than  that  of  simultaneously  observed  Evershed 
flows.  vzc  and  V[p  are  generally  well  correlated  with  z/,-.  The 
linear  fittings  show  that  \vip\  and  zz,  are  close  to  each  other 
and  much  smaller  than  \vzc\.  The  large  \vzc\  may  be  ar¬ 
tificial  signal  resulted  from  large  Stokes  V  asymmetry  and 
limited  spectral  resolution  as  discussed  above.  It  is  also 
possible  that  the  mixed  polarities  there  contain  a  compo- 
nent  possessing  a  high  speed  downflow,  as  observed  earlier 
by  Martinez  Pillet  et  akl  (119941).  Detailed  examination  of  the 
large  redshift  requires  two-component  Stokes  inversion  of  in¬ 
dividual  profiles,  which  is  out  of  the  scope  of  this  study. 

4.3.  Asymmetry  of  Stokes  Profiles 

Asymmetries  of  both  Stokes  V  and  LP  profiles  were  cal¬ 
culated.  We  present  results  from  Stokes  V  only  because  the 
noisier  and  weaker  LP  signals  produce  large  uncertainty  in  de¬ 
termined  asymmetries  although  their  spatial  distributions  are 
qualitatively  similar  to  those  of  Stokes  V.  Figure  [5]  shows  the 
spatial  distribution  and  spectral  variation  of  the  Stokes  V  am¬ 
plitude  asymmetry  (8a)  and  area  asymmetry  (<5A).  We  also 
plot  their  histograms  for  different  magnetic  regions  in  Fig¬ 
ure  [6]  in  order  to  quantitatively  compare  the  Stokes  V  asym¬ 
metries  among  the  three  spectral  lines.  The  typical  Stokes  I,V 
and  LP  profiles  in  those  regions  are  also  plotted  in  Figure [7] 
We  again  discuss  various  magnetic  features  separately  as  fol¬ 
lows. 

For  plages,  the  red  lobes  are  broader  and  shallower  than 
the  blue  lobes  for  the  two  photospheric  Fe  I  lines,  which 
makes  highly  asymmetric  Stokes  V  profiles  thus  substan¬ 
tial  positive  8a  and  8A  with  8a  >  8A  generally,  as  pre¬ 
viously  found  (IStenflo  et  all  1 19841  ISolanki  &  Stenflol  11984 
1 19851  iGrossmann-Doerth  et  alJl!996h  ISigwarthll200H).  It  is 
believed  that  the  “magnetic  canopy”  effect  and  the  small 
MFF  in  photospheric  plages  play  important  roles  in  pro¬ 
ducing  such  large  positive  asymmetries  (e.g.jLeka  &  Steinerl 
120011).  A  noteworthy  finding  here  is  that  the  chromospheric 
Mg  I  Stokes  V  profiles  are  much  more  symmetric  as  can  be 
seen  that  the  most  probable  8a  and  8A  are  very  small  posi¬ 
tive  values.  iBriand  &  Solankil  (119981)  also  observed  smaller 
Stokes  V  asymmetry  in  the  Mg  I  bo  line  compared  to  a  photo¬ 
spheric  line  in  network  regions.  It  was  found  that  for  small  or 
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Fig.  4. —  Upper  panels :  Dopplergrams  made  using  Stokes  /  line  center  shifts  (i/j),  Stokes  V  ZC  shifts  (vzc)  and  LP  profiles  shifts  (vip)  for  the  three  spectral 
lines.  The  numbers  in  unit  of  km  s-1  are  the  mean  velocities  within  the  strong  and  compact  redshift  region  around  the  MPIL  of  the  5  spot,  and  those  in  the 
penumbral  region  outlined  by  the  white  box.  The  spatial  points  with  noisy  and  abnormal  profiles  are  assigned  a  value  of  0.  Note  that  the  methods  of  determining 
Vi,  vzc  and  v/p  may  break  down  in  some  umbral  areas  for  fully  split  lines.  Lower  panels :  Scatter  plots  of  Vj  vs.  vzc  (forth  row )  and  v,  vs.  Vjp  (fifth  row)  for 
each  spectral  line.  The  color-coded  data  points  stand  for  different  magnetic  regions.  The  slope  k  and  the  linear  correlation  coefficient  (C.C.)  from  the  line  fits 
vzc  =  kvi  +  b  or  V[p  =  kvj  +  b  for  each  group  of  data  are  labeled. 
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Fig.  5. —  Images  constructed  with  values  of  Stokes  V  amplitude  asymmetries  {5a,  first  row)  and  area  asymmetries  ( 5A ,  second  row)  at  each  spatial  point  for 
the  three  spectral  lines.  The  magnitude  of  asymmetry  is  indicated  by  the  scales  at  right.  The  spatial  points  with  noisy  and  abnormal  profiles  are  assigned  a  value 
of  0. 
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Fig.  6. —  The  distribution  of  Stokes  V  amplitude  asymmetries  (5a,  first  row)  and  area  asymmetries  (5A,  second  row)  in  different  magnetic  regions  for  the  three 
spectral  lines.  The  black,  gray,  and  dotted  lines  are  for  the  Fe  I  630.25,  Fe  1  630.15,  and  Mg  1 517.27  observations,  respectively. 


weak  magnetic  structures,  such  as  plage  and  network,  the  de¬ 
gree  of  Stokes  asymmetries  decreases  if  MFF  increases  (e.g., 
ZaveretaT]|1990t  iMartlnez  Pillet  et  all  1 19971  ISigwarth  et  all 
1999f).  The  dependence  of  amplitude  of  Stokes  V  asymme¬ 
tries  on  some  line  parameters,  such  as  Doppler  shift.  Zee- 
man  shift,  line  strength  and  width,  has  also  been  heuristically 
considered  within  the  frame  of  a  simplified  magnetic  canopy 
model  dGrossmann-Doerth  et  al.lf]~989l).  There  the  asymme¬ 
try  is  expected  to  be  small  for  a  same  downdraft  velocity  if 
the  line  is  broad  and  not  strongly  Zeeman  split.  Our  result 
confirms  these  trends  from  observational  point  of  view.  To 
sum  up,  both  the  broader  line  width  and  the  larger  MFF  in  the 
formation  layer  of  the  Mg  1 1>2  line  are  responsible  for  the  ob¬ 
served  small  asymmetry.  Together  with  the  aforementioned 
finding  that  some  of  the  conspicuous  downflows  surrounding 
the  plages  disappear  in  the  Mg  I  v,  Dopplergram,  we  there¬ 
fore  suspect  that  the  “magnetic  canopy”  of  plages  could  be  a 
shallow  phenomenon  mainly  below  the  TMR,  which  however 
needs  more  solid  observational  evidence  to  justify. 


For  penumbrae  (especially  their  outer  parts)  and  pores,  5a 
and  5A  are  primarily  positive  and  negative  in  the  photosphere 
and  in  the  low  chromosphere,  respectively.  Previous  studies 
have  shown  that  penumbrae  and  pores  share  some  common 
properties.  They  both  present  an  interface  with  the  quiet  sun 
and  are  surrounded  by  vortex  convective  flows  with  similar 
patterns  (e.g.lDeng  et  alJl2007l).  This  study  hence  adds  to  the 
evidence  of  another  common  property,  i.e.,  Stokes  V  asymme¬ 
tries  and  their  variation  with  height.  Unlike  photospheric  lines 
that  can  be  easily  synthesized  under  LTE  conditions,  the  phys¬ 
ical  origins  of  the  asymmetries  of  chromospheric  lines  formed 
in  Non-LTE  atmosphere  are  far  from  fully  studied.  Here  we 
can  only  present  a  speculation  for  the  sign  reversal  of  Stokes  V 
asymmetries  based  on  a  logical  analog  of  chromospheric  line 
asymmetries  to  the  better  known  photospheric  line  asymme¬ 
tries  that  are  mainly  caused  by  velocity  gradients.  The  change 
in  sign  of  Evershed  effect  for  penumbrae  may  be  responsible 
for  some  but  not  all  areas,  because  the  Stokes  V  asymme¬ 
tries  reverse  sign  all  over  the  outer  penumbral  borders  while 
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the  inverse  of  Evershed  flow  only  occurs  in  some  penumbral 
sectors.  We  speculate  that  the  sign  reversal  of  Stokes  V  asym¬ 
metries  in  both  outer  penumbrae  and  pores  could  indicate  that 
the  velocity  in  and  around  them  may  change  drastically  from 
the  photosphere  to  the  low  chromosphere.  Figure 0] provides 
clues  to  some  of  the  changes  because  some  downflows  in  and 
around  pores  do  not  hold  in  the  Mg  I Dopplergram  instead 
more  upflows  can  be  seen  (e.g.,  in  the  scatter  plot  of  the  Mg  I 
line).  Moreover.lSankarasubramanian  &  Rimmelel  d2003l)  ob¬ 
served  pores  in  different  wavelengths  and  found  that  veloci¬ 
ties  in  the  narrow  downflow  region  around  pores  tend  to  de¬ 
crease  with  height  at  first,  and  then  change  to  upflow  in  the 
low  chromosphere.  Such  tendency  is  also  visible  in  the  MHD 
simulation  of  pores  (see  Fig.  3  oflLeka  &  Steineril200ll).  The 
Stokes  V  asymmetries  change  from  positive  for  the  Fe  I  lines 
to  negative  for  the  Mg  I  line,  if  indeed  all  line  asymmetries 
have  a  velocity  gradient  origin,  thus  provides  important  clues 
on  the  3D  structure  of  the  flow  field  in  and  around  pores  and 
outer  penumbrae. 

For  umbrae,  since  they  lack  of  mass  flows  and  have  MFF  of 
~1  in  the  low  photosphere,  both  of  the  most  probable  5a  and 
5A  of  the  Fe  1 630.25  nm  line  are  close  to  0.  The  large  Zeeman 
splitting  of  the  line  in  strong  field  may  also  contribute  to  the 
observed  small  asymmetry  (iGrossmann-Doerth et  al.|il989). 
However,  the  distributions  of  5a  and  5A  of  the  other  two  spec¬ 
tral  lines  show  a  clear  shift  toward  negative.  Besides  the  noise 
induced  inaccuracy  of  the  determined  asymmetry  due  to  low 
light  level  in  the  umbra,  it  is  also  possible  that  as  height  in¬ 
creases  umbrae  become  less  uniform  and  static  compared  with 
their  structure  in  the  low  photosphere. 

For  the  compact  region  with  strong  redshifts  around  the 
MPIL  of  the  5  spot,  both  5a  and  5A  have  prominent  values  that 
increase  with  height.  It  seems  from  Figure[7]that  the  red  wings 
of  the  Stokes  /,  V,  and  LP  profiles  of  the  two  photospheric 
lines  obviously  have  further  redward  extension.  Meanwhile, 
the  blue  damping  wings  of  the  Mg  I  517.27  nm  line  are  sur¬ 
prisingly  enhanced,  which  remains  a  puzzling  issue.  We  con¬ 
jecture  that  the  observed  large  asymmetries  could  be  due  to 
the  combination  of  the  following  two  effects:  the  mixed  po¬ 
larities  or  flows,  and  the  large  gradients  in  velocity,  magnetic 
field  vector,  or  temperature  near  the  MPIL  region.  Consider¬ 
ing  all  the  unusual  properties  observed  in  this  region,  we  sug¬ 
gest  that  the  magnetic  and  flow  fields  around  the  MPIL  of  the 
5  spot  are  very  complex,  which  deserve  a  deeper  investigation 
in  the  future. 

5.  SUMMARY 

We  have  presented  a  comprehensive  study  of  the  Doppler 
shifts  and  asymmetric  properties  of  the  Stokes  profiles,  the 
uniqueness  of  which  is  especially  in  the  combination  of  mag¬ 
netic  sensitive  spectral  lines  that  form  at  different  heights 
spanning  from  the  photosphere  to  the  low  chromosphere. 
Our  quantitative  analysis  benefits  from  the  simultaneous  and 
high  resolution  spectro-polarimetric  observation  that  include 
a  chromospheric  line.  We  summarize  the  major  findings  and 
interpretations  as  follows. 

1 .  Vi  is  very  close  to  v/p  but  significantly  different  from  vzc 
in  penumbral  areas  near  disc  center  for  the  photospheric 
Fe  I  lines,  which  provides  direct  and  strong  evidence 
that  the  penumbral  Evershed  flows  are  magnetized  and 
mainly  carried  by  the  horizontal  magnetic  component. 

2.  The  rudimentary  inverse  Evershed  effect  observed  by 


the  Mg  I  b2  line  provides  a  qualitative  evidence  on  its 
formation  height  that  is  around  or  just  above  the  TMR. 

3.  vzc  and  vip  in  penumbrae  and  vzc  in  pores  generally  ap¬ 
proach  their  observed  by  the  chromospheric  Mg  I 
line,  which  is  not  the  case  for  the  photospheric  Fe  I 
lines.  We  speculate  that  the  interlocking-comb  struc¬ 
ture  of  the  penumbral  magnetic  and  flow  field  might  be 
a  shallow  phenomenon  mainly  in  the  photosphere. 

4.  The  Stokes  V  asymmetries  exhibit  similar  behavior  in 
the  regions  of  outer  penumbrae  and  pores.  Both  5a  and 
5A  evolve  from  primarily  positive  values  observed  by 
photospheric  Fe  I  lines  to  primarily  negative  values  ob¬ 
served  by  the  chromospheric  Mg  I  line. 

5.  In  plage  regions  where  the  “magnetic  canopy”  effect 
is  present,  Stokes  V  profiles  are  highly  asymmetric  in 
the  photosphere  but  more  symmetric  in  the  low  chro¬ 
mosphere.  Moreover,  some  of  the  photospheric  down¬ 
flows  surrounding  the  plages  disappear  in  the  Mg  I 
Dopplergram.  These  suggest  that  the  plage  “magnetic 
canopy”  could  be  a  shallow  phenomenon  mainly  below 
the  TMR. 

6.  All  the  Vi,  v-c  and  V[p  Dopplergrams  reveal  strong  red- 
shifts  in  a  compact  region  around  the  MPIL  and  within 
the  common  penumbra  of  the  5  spot.  Large  posi¬ 
tive  Stokes  V  asymmetries  are  also  present  and  even 
strengthen  with  height.  These  unusual  phenomena  are 
manifestations  of  the  complex  nature  of  the  3D  topol¬ 
ogy  of  the  highly  non-potential  region. 

We  conclude  that  analysis  of  Stokes  profiles  using  both 
the  photospheric  and  chromospheric  spectral  lines  can  fur¬ 
nish  crucial  information  on  the  3D  structure  of  the  magnetic 
and  flow  fields.  Similar  study  using  other  chromospheric  lines 
will  be  essential  to  see  if  the  behavior  and  trend  of  velocities 
and  asymmetries  observed  by  the  Mg  I  £>2  line  is  general  for  all 
the  chromospheric  lines.  The  speculations  made  in  this  study, 
such  as  the  velocity  gradient  origin  of  the  chromospheric  line 
asymmetries  and  the  thinness  of  the  interlocking-comb  struc¬ 
ture  of  penumbral  magnetic  and  flow  fields,  need  further  in¬ 
vestigation  to  justify.  To  advance  firmly  towards  a  sound  and 
reliable  measurement  of  magnetic  fields  in  the  chromosphere, 
not  only  more  polarimetric  observations  using  chromospheric 
lines  are  needed,  but  also  a  better  understanding  of  their  com¬ 
plicated  profile  formation  conditions  is  critical. 
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